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Abstract 

We have computed the full next-to-leading (NLO) QCD corrections to the differential 
distributions d 2 a/ (dpx dy) for pseudo-scalar Higgs (A) production at large hadron colliders. 
This calculation has been carried out using the effective Lagrangian approach which is valid 
as long as the mass of the pseudo-scalar Higgs boson tba and its transverse momentum px do 
not exceed the top-quark mass rat- The shape of the distributions hardly differ from those 
obtained for scalar Higgs (H) production because, apart from the overall coupling constant 
and mass, there are only small differences between the partonic differential distributions for 
scalar and pseudo-scalar production. Therefore there are only differences in the magnitudes 
of the hadronic differential distributions which can be mainly attributed to the unknown 
mixing angle (3 describing the pseudo-scalar Higgs coupling to the top quarks. 
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The scalar Higgs boson H, which is the corner stone of the standard model, is the only particle 
which has not yet been observed. Its discovery or its absence will shed light on the mechanism 
how particles acquire mass as well as answer questions about super-symmetric extensions of the 
standard model or about the compositeness of the existing particles and the Higgs boson. Among 
these two alternatives super-symmetry is the most appealing one, in particular the minimal 
super-symmetric extension of the standard model. The latter version contains two complex 
Higgs doublets instead of one and it is therefore called the Two-Higgs-Doublet Model (2HDM). 
Here the scalar particle spectrum contains both the Higgs boson H and another neutral scalar 
boson h. Furthermore it contains two charged scalar bosons H 1 * 1 and a neutral pseudo-scalar 
Higgs boson A. The tree-level masses are expressed in two independent parameters, namely the 
mass wa and the ratio of the vacuum expectation values of the two Higgs doublets defined by 
tan/? = V2/V1 (see e.g. Jl|]). According to the experiments at LEP their parameter ranges are 
restricted so that m A < 91.9 GeV/c 2 and 0.5 < tan/3 < 2.4 || are excluded. In this paper we 
study A-production which in lowest order proceeds via gluon-gluon fusion where the gluons are 
coupled to the A via a heavy flavour triangular loop. This is similar to H-production except 
that now the coupling constant describing the interaction of the A with the quarks depends on 
both the masses of the quarks and on the angle (3. This follows from the 2HDM where the 
coupling constants of the up and down quarks behave like g up ~ m u cot (3 and gdown ~ m d tan (3 
respectively |IJ • Since the effective Lagrangian approach below is only valid in the case the mass 
of the quark appearing in the triangular loop satisfies the condition m q ^> mA, the bottom quark 
is excluded. However then we have to require that in the 2HDM the coupling of the A to the 
top-quark is stronger than to the bottom-quark which implies the condition 

— >tan 2 /3. (1) 
m b 

If we choose m b = 4.5 GeV/c 2 and m t = 173.4 GeV/c 2 one obtains the inequality tan/5 <C 6.21. 
In view of the experimental boundaries above one can conclude that the results of the calculation 
below can be only applied for the regions tan/3 < 0.5 and 2.4 < tan/3 < 6.21. 

In the effective Lagrangian approach scalar H-production is described by the Lagrangian 
density §, 

C? ff = G n <S> R (x)0(x), with 0(x) = -~G« u (x)G a ^(x), (2) 
whereas pseudo-scalar A-production is obtained from 0, ||, 



G A 1 (x) + G A 2 (x) 



with 
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O^x) = -- e, uXa G?(x) Gl°{x) , Q 2 (x) = £ qi{x) W q^x) , (3) 



=1 



where $ H (x) and $ A (x) are the scalar and pseudo-scalar fields respectively and rif denotes the 
number of light flavours. Up to NLO the operator O^x) only contributes when it interferes with 
the operator Oi(x) provided the quarks are massless. The effective couplings Gb (B = H, A) are 
determined by the top-quark triangular graph describing the decay process B — > g + g in the 
limit m t — ► 00 

G| = 4v^f^VG F r 2 F 2 (r B )C 2 L( /U 2 ),4^ r B = , B = H,A. (4) 
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and the functions Fb are defined by 

F H (r) = 1 + (1 - r) /(r) , F A (r) = /(r) cot/? 

fir) = arcsin 2 —= , for r > 1 , 




/(r) = --lln - ; v ^ - + 7tz) , for r < 1 . (V) 

In the large mj- limit F(r) behaves as 

lim F H (r) = ^- , lim F A (r) = - cot /3 . (6) 

t— >oo O T r^oo T 

Here m and m t denote the masses of the (pseudo-) scalar Higgs boson and the top quark respec- 
tively. The running coupling constant is given by a s (n 2 ) where \i r denotes the renormalization 
scale and Gp is the Fermi constant. The coefficient functions Cb originate from the corrections 
to the top-quark triangular graph provided one takes the limit m t — > oo. We have presented the 
couplings Gb in Eq. (4) for general m t on the Born level only in order to keep some part of the 
top-quark mass dependence. This is an approximation because the gluons which couple to the 
(pseudo-) scalar Higgs boson via the top-quark loop in the partonic subprocesses are very often 
virtual. The virtual-gluon momentum dependence is neither described by Fb(t) nor by C-q- For 
on-mass-shell gluons the latter quantity has been computed in the large mt limit up to order 
a s in ||, [|], || and up to order a 2 in j4j, [7). For our NLO calculations we only need these 
coefficient functions corrected up to order a s and they read 

Cu(a s (^ r ),^) = l + + (7) 



C A a s (i4 
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where is presented in a five-flavour number scheme. Notice that Eq. (8) holds in all 
orders because of the Adler-Bardeen theorem ||. The effective Lagrangian approach has been 
successfully applied to compute the total cross section of scalar Higgs production in hadron- 
hadron collisions in NLO and NNLO @, 0, Q, @, 0. In the case of pseudo-scalar 



Higgs production this cross section was computed in NLO in H and in NNLO in [Tl|, [P^. 

In this paper we study the semi-inclusive reaction with one pseudo-scalar Higgs boson A in 
the final state which is given by 

HifrO + H 2 (p 2 ) - A(-p 5 ) +' X' , (9) 

where Hx and H 2 denote the incoming hadrons and X represents an inclusive hadronic final 
state. Further we define the following kinematical invariants 

S=( Pl + P2 ) 2 , T=( Pl+P5 ) 2 , u=( P2 +p 5 ) 2 . (10) 

The latter two invariants can be expressed in terms of the transverse momentum px and rapidity 
y variables as 



T = m 2 — \/Pt + m 2 cosh y + v^S \jp\ + m2 s i nn V 



U = m 2 — ^/S J p\ + m 2 cosh y — \jp\ + m 2 sinh y , (11 
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where m is the mass of the pseudo-scalar Higgs boson. The hadronic cross section is given by 



d 2 rr HlH2 



X\ Jx 2 ,min x 2 at a u 

—U — X\(T — m 2 ) — m 2 . s 

with xi >min = — — , x 2Mn = — — — , (12) 

b + 1 — mr X\b + U — mr 

where s, t and u are the partonic analogues of S, T and U in Eq. (10) where pi and p 2 now 
represent the incoming parton momenta. Further f^' 1 denotes the parton density corresponding 
to hadron H; and /i stands for the factorization scale which for convenience is set equal to the 
renormalization scale fi r appearing in Eq. (4). The NLO corrections to the partonic cross 
section d 2 a/(dt du) in the case of H-production based on the effective Lagrangian in Eq. (2) are 
presented in [[U| and |ji"7| . Here we will give the corresponding results for the A described by 
the Lagrangian in Eq. (3). The calculation proceeds in the same way as presented in flL7|| . We 
use n-dimensional regularization in order to compute the loop and phase space integrals which 
contain ultraviolet, infrared and collinear singularities. However there is one extra complication 
in the pseudo-scalar case. This concerns the Levi-Civita tensor in Eq. (3) which is essentially a 



four dimensional object. Here we follow the same prescription as in Eq. (4) in |L4|, |L5| where 
the product of two Levi-Civita tensors is contracted in n dimensions if one sums over dummy 
Lorentz indices. This prescription leads to an interference between diagrams carrying the vertex 
coming from the operator 0\{x) and those carrying the vertex corresponding to the operator 



2 (x) in Eq. (3) (see 0]). The LO subprocesses contributing to the partonic cross section are 
given by 

g + g^g + A, q + q^g + A, q{q) + g — q{q) + A . (13) 

The matrix elements squared do not differ from those derived for the scalar H provided n = 4, 
see 0, 0, which implies that the LO double differential partonic cross sections are the same 
for both bosons except for an overall constant given by F^(t) in Eq. (5). In NLO one has to 
compute the one-loop virtual corrections to the processes in Eq. (13) above and to add the 
contributions from the following two-to-three-body reactions 

g + g^ g + g + A, g + g^q. + qi + A, (14) 

Q + Q 9 + 9 + A, qi + q 2 -> <?i + q 2 + A q 1 ^ q 2 , 

q + q^qi + qi + A q, t ^ q , q + q -> q + q + A, (15) 

Ql + Q2 -> qi + Q2 + A q!^q 2 , q + q q + q + A , (16) 

q(q)+g ^q(q)+g + A. (17) 

After renormalization of the strong coupling constant a s and mass factorization which are car- 
ried out in the MS-scheme we obtain the NLO corrected coefficient functions according to the 



procedure in ||17|| . The coefficient functions are as long as in the case of H-production so that 
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they cannot be explicitly presented. However the differences between the results for the H and 
the A are so small that we can show them below. If we put for simplicity Gn = Ga = G and 
mn = Ota = m the differences between the soft-plus-virtual differential cross sections are given 
by 



d 2 S+V j2 s+v 

c 2 a a gg^g A g2 Q gg^g H 



d td u 



d td u 



nS(s + t + u-m 2 )G 2 (^4^ 
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gg^g Bl 



j2 ^s+v 7 2 s+v 

g 2 " U qq^g A U qq-*gU 



d td u 



d td u 



2CJM 



(i) 



qq^g Bl 



Cw-Ca 1MB, 



(i) I 

99-»9 Bl 



(19) 



J2 S+V -2 ^.S+V 

^2 " u gg—*g A g 2 u qg^q H 



d td u 



d td u 



n5(s + t + u-m 2 )G 2 



47T 



iV(iV 2 



2C A b 



Cw-Ca WMB 



(i) I 

99^9 Bl 



(20) 



where we have added to the righthand side of Eqs. (19) and (20) the contributions coming from 



the interference of the graphs in Figs, lb, Id with those in Figs. 2b, 2d which are shown in |fL4 
The colour factors of the group SU(N) are given by Ca = N and Cp = (N 2 — 1)/(2N) and the 
Born matrix elements squared belonging to the processes in Eq. (13) are equal to 



IM (1) I 2 

\ 1V1 gg^g Bl 



iV(iV 2 - 1) 



stu 



s A + t 4 + m 4 + m* 



(21) 



IM (1) I 2 

\ m qq->g Bl 



Ca C f 



t 2 + u 2 



(22) 



IM (1) I 2 



Ca Cf 



u 



-s 2 -t 2 



(23) 



The differences above can be wholly attributed to the virtual corrections and not to the soft 
gluon contributions which are the same for both H and A-production. These virtual corrections 
also entail some extra terms denoted by 



1MB 
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gg^g Bl 



-N(N 2 
3 V 
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stu 



stu + m 2 (st + su + t 
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(24) 
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\MB^ gB \ 2 = C A C F (-t-u), 



1MB 



(i) I 

qg^q B\ 



= C A C F (s + t). 



Denoting the two-to-three-body reactions by 

a(pi) + b(p 2 ) -> c(-p 3 ) + d(-p 4 ) + A(-p5) , s 4 = (p 3 + p 4 ) 2 , 

then the differences between the partonic cross sections due to the subprocesses in Eqs. 
are equal to 

j2 HARD j2 HARD 

„2 U 99^99 A _ 2 " 99^99 H _ 



(25) 
(26) 

(27) 
(14)-(17) 
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(36) 



where the meaning of the superscript HARD is explained in |[L7| . From these expressions we 
infer that the partonic cross sections (coefficient functions) for H and A are equal in LO and 
almost equal in NLO. This means that apart from the overall normalisation due to the constant 
Gb there will not be any difference in the shapes of the double differential cross sections. We 
show this in Fig. 1 where we plot the ratio 



R 



d&A 
da-R 



(37) 



for d<j g = dus/dpr and proton-proton collisions at the LHC with \/~S = 14 TeV. For these 
and the next plots we have adopted the parton density set MRST98 (LO, lo05a.dat) [|l8j for 
the LO calculations with A^ LO = 130.5 MeV as input for the leading log running coupling 
constant. For the NLO cross sections we have chosen the set MTST99 (NLO, cor01.dat) ]19|] with 
A^ LO = 220 MeV as input for the next-to- leading log running coupling constant. Furthermore 
the factorization/renormalization scale is chosen to be fi 2 = fi 2 = p 2 ^ + Mg. For the masses 
of the Higgs bosons we take mn = rriA = 120 GeV/c 2 and the top quark mass is set equal to 
m t = 173.4 GeV/c 2 . Further we have put tan/? = 1. In the case of an infinite top quark mass, 
(here we choose m t = 173.4 x 10 3 GeV/c 2 ), we get R h ° = 9/4 irrespective of the values of 
mn and m^. This follows from Eq. (6) and the fact that the LO partonic cross sections are 
the same for H-production and A-production. A finite m t as given above introduces a small 
effect and one gets R LO = 2.31 which amounts to a shift upwards of 0.06 (see Fig. 1). In NLO 
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the partonic cross sections differ a little bit and C\ = [1 + 22a s / (47r)] C\ (see Eqs. (7), (8)). 
Therefore we expect a deviation from the R LO result when m t is taken infinite in both the LO 
and NLO reactions. However it turns out that both differences compensate each other. The 
NLO corrected partonic cross section for A is larger than the one for H and one obtains an 
upward shift Ai? NLO = 0.26. The shift due to the coefficient function in Eq. (7) is negative and 
amounts to Ai? NLO = —0.24. Hence the actual value becomes i? NLO = 2.27 (see Fig. 1) which 
is very close to R NLO = 9/4. If is finite one gets again an upward shift of 0.06 like in LO 
and i? NLO = 2.33 (see Fig. 1). One can make similar plots for the rapidity y distributions which 
yield the same ratios as shown in Fig. 1 for the pr distributions. The most important feature is 
that the ratios are independent of pt and y showing the shape independence of the distributions 
on the parity of the Higgs boson (scalar versus pseudo-scalar). This behaviour was discovered 
for both the (pseudo)-scalar p? distributions and for the opening angle distribution between the 
(pseudo)-scalar bosons and the highest py-jet in the reaction p + p — * (H or A) + jet + jet +' X' 
in [pOfl . From Fig. 1 and the observations made above it is clear that the ratios between 
the NLO and LO corrected cross sections (K-factors) are the same for H-production and A- 
production. This also holds for the variation of the NLO cross sections with respect to the mass 
factorization/renormalization scales. They are given for H-production in [[17] and we do not 
have to show them again for A-production. In Fig. 2 we present the pt distributions in NLO for 
A-production in proton-antiproton collisions at \/S = 2 TeV (Fermilab Tevatron, Run II) and 
in proton-proton collisions at v^S* = 14 TeV (LHC). Further we have chosen m& = 91.9 GeV/c 2 
and tan /3 = 0.5. The parton density set and the factorization scale are given above. From Fig. 2 
we infer that the p^-distributions decrease rather slowly as pt increases and that the differential 
cross section for the Tevatron is two orders of magnitude smaller than the one predicted for the 
LHC. The latter observation also holds for the corresponding rapidity distributions shown in 
Fig. 3. They are obtained by integrating d 2 a^/{dp T dy) over the range PT,mm < Pt < 8 Pt, mm 
with pr.min = 30 GeV/c. The cross section for px > 8 PT,mm is negligible. Notice that the range 
of the rapidity for A-production at the Tevatron is rather small. Finally we want to comment 
on the relative importance of the partonic subprocesses contributing to the hadronic differential 
cross section in Eq. (12). For the LHC (y/S = 14 TeV) the gg-channel dominates and the 
gg-subprocess contributes about one third of the cross section. This is because at these high 
energies the x- values of the gluon density fg(x) is so small that it becomes much larger than the 
quark densities. At lower energies like \/S = 2 TeV (Tevatron) the x-values are larger so that the 
valence quark densities also play a role. This explains why the contribution of the gg-subprocess 
is of the same magnitude as the one from the gg-channel for A-production at the Tevatron. 
Acknowledgments. We are indebted to V. Ravindran for discussions about the relevance of 
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Figure Captions 
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Fig. 1. The ratio R in Eq. (37) plotted as a function of pt for v^S* = 14 TeV and /i 2 = p\ + 

with m n = m A = 120 GeV/c 2 ; B LO (m t = oo) (dotted line), R LO (m t = 173.4) (solid line) 
R NLO (m t = oo) (dot-dashed line) R NLO (m t = 173.4) (dashed line). 

Fig. 2. The transverse momentum distribution da^/dpr with fi 2 = p T +m\, uia = 91.9 GeV/c 2 , 
tan/3 = 0.5; VS = 14 TeV (solid line), y/S = 2 TeV (dashed line). 

Fig. 3. The rapidity distribution da^/dy calculated from the integral of d 2 a^/ (dpr dy) between 
8 Pr,min > Pt > PT,min and pr,min = 30 GeV/c. Input parameters are /x 2 = p Tra ^ + m\, 
m A = 91.9 GeV/c 2 , tan/5 = 0.5; VS = 14 TeV (solid line), VS = 2 TeV (dashed line). 
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